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Abstract. Transverse mass spectra of pions and protons measured in central collisions of heavy ions at the SPS 
and at RHIC are compared to a hydrodynamic parameterization. While the chemical temperature needed at RHIC 
is significantly higher compared to SPS, the spectra may be described using kinetic freeze-out parameters which are 
similar for both beam energies. At RHIC either the temperature or the flow velocity is higher, but the data provide 
no unambiguous proof for much stronger transverse flow. The contribution of such hydrodynamic emission at high 
transverse momenta is investigated in detail. It is shown that hydrodynamics may be relevant up to relatively high 
transverse momenta. The importance of the velocity profile used in this context is highlighted. 

PACS. 25.75.Dw Particle and resonance production 



1 Introduction 

The description of transverse momentum (or transverse mass) 
spectra with hydrodynamic models is one of the well estab 
lished tools in ultrarelativistic heavy-ion physics 
The applicability of such models is still under discussion, be- 
cause they require local thermalization, which may not hold 
for the entire reaction volume. Furthermore, at high transverse 
momenta hard scattering is expected to provide an important 
contribution to particle production - the exact momentum re- 
gion where this would apply is under debate. 

Comparisons to experimental data have often been per- 
formed using parameterizations of hydrodynamic distributions, 
which have to use simple assumptions about the freeze-out 
conditions, e.g. mostly one universal freeze-out temperature is 
used. While such parameterizations can not replace full hydro- 
dynamic calculations, they can provide a reasonable guideline 
to the general behavior of particle distributions, especially as 
systematic studies are more accessible in this case due to the 
much more modest calculational effort. 

In the present paper I will assume the validity of hydrody- 
namic pictures, but I will not attempt to deliver a full-fledged, 
reaUstic hydrodynamic calculation. A parameterization will be 
used to describe momentum spectra at low and intermediate 
rriT- f\ Parameters for chemical and kinetic freeze-out will be 
determined from a simultaneous fit of pion and (anti-)proton 
spectra. The main intention is to use a not too unrealistic hy- 
drodynamic parameterization fitted to low transverse masses to 
study the extrapolation to high rriT- 

This is very important also in the light of recent attempts to 
measure jet quenching at RHIC via transverse momentum dis- 



tributions [^. If the particle emission at low mr is governed 
by hydrodynamics, as e.g. the results on elliptic flow [ p^pl] ] 
suggest, then such a hydrodynamic source will inevitably con- 
tribute at high rriT also. Such a hydrodynamic contribution 
would have to be taken into account before comparing them 
to pQCD calculations and obtaining estimates of the amount of 
suppression of hard scattering. 



2 The Model 

The calculations presented in this paper are based on a parame- 
terization of the source at freeze-out by Wiedemann and Heinz 
iQl which is motivated by a hydrodynamic approach and in- 
cludes effects of transverse flow and resonance decays. The 
original computer program calculates the direct production of 
pions and the contributions from the most important resonances 
having two- or three-body decays including pions (p, K^, K*, 
A, S + A, 7y, u!, rj'). The transverse momentum spectra of the 
directly emitted resonances r are given by: 

j2" = const. • {2Jr + 1) 



dM4 



Jo 



(1) 



' While such calculations can also be used to compare to two- 
particle-correlations, this will not be the scope of this paper 



where ^ — r/ R and R is the Gaussian radius of the source. Be- 
cause the integrals are evaluated numerically, fixed integration 
boundaries have to be chosen. The upper limit of ^ = 4 was 
chosen to ensure that for most reasonable shapes only negligi- 
ble tails of the distribution extend beyond the limits. The trans- 
verse rapidity of the source element is given as 77t(C) = ?7/C"' 
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where the default value of the power is n = 1. The transverse 
rapidity parameter 77 f controls the amount of transverse flow. 

This distribution represents the rapidity integrated spectrum, 
which will be used in this paper, as it is easier to calculate 
compared to the general rapidity differential spectrum. In com- 
parison to experimentally measured spectra, which are mostly 
from limited rapidity regions, this may introduce a bias. This 
bias should however be small, if the rapidity dependence of 
the spectra is negligible as in a longitudinal scaling expansion 
(Bjorken) scenario. In fact, although the Bjorken picture is not 
at all applicable at SPS, the dependence of the spectral shape 
on rapidity is already reasonably small [|l2||. The dependence 
of yields of different species on rapidity may have a small ef- 
fect on the chemical parameters in the fits. We will comment 
on this again in a later section. 

The original version of this model uses the foUowing as- 
sumptions: 

1 . The spatial density distribution at freeze-out is chosen as a 
Gaussian. 

2. One universal freeze-out temperature is used which deter- 
mines both the spectral shape and the ratio of different par- 
ticle species (i.e. resonances). 

In the calculations performed here the model has been mod- 
ified in the following way: 

1. As the spatial distribution a Woods-Saxon shape: 



Pws(0 = PQ ■ 



1 



l + cxpZ\(e-l) 



(2) 



has been introduced. By varying the parameter A the shape 
can be adjusted. For large values of A the shape approaches 
a box-like distribution. The shape may also be chosen very 
similar to a Gaussian for comparison. The spatial distribu- 
tions fix implicitly the shape of the velocity distribution. 
While this may not be very relevant at low transverse mo- 
menta, it has l3een shown to be important at high px [^. 
The influence of different distributions has also been inves- 
tigated in |§. 

2. The spectra for a given particle (or resonance) are calcu- 
lated using the kinetic temperature Tkin- Finally the nor- 
malization is readjusted to the chemical temperature Tchem 
assuming that dN/dy at midrapidity scales with the tem- 
perature as nifl: 



dy ~ (27r)2^ KT"^ 



V ^fm^ 



2™ 



2 cxp - 



T 



(3) 



It is commonly believed that chemical freeze-out (deter- 
mining particle ratios) should occur at higher temperature 
than kinetic freeze-out (determining the spectral shape), so 
at least two independent temperatures may be needed. 
3. Furthermore the program has been enhanced to simultane- 
ously describe protons and antiprotons in addition to pions. 
Here the decay contributions from A and S + A have been 
taken into account. This of course requires the introduction 
of another parameter: the baryonic chemical potential jis- 

To better understand the importance of the different shapes 
of the source distribution, it is helpful to investigate the result- 
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Fig. 1. Velocity distributions of the model source for different source 
shapes at an average velocity of {Pt) ~ 0.12. 



ing velocity distributions: 



dN 
dp^ 

and their mean values: 

(/3t) 



1- [Prior 



(4) 



(5) 



as they directly influence the shape of the momentum distribu- 
tion, while the spatial distribution is actually not directly rele- 
vant. Fig. [l| shows velocity distributions for different shapes at 
a small average velocity of {(3t) ~ 0.12. For such small veloc- 
ities the relation between rapidity and velocity fir = tanh % 
may be approximated by a linear relation, so most of the ex- 
pected properties of the distributions are preserved. One can 
see that the box profile yields a sharp cutoff of the velocity dis- 
tribution. The Woods-Saxon with Z\ = 50 looks very similar 
to the box with a slightly smeared out edge. The Gaussian has 
a much more smooth edge resulting from the tail of the spatial 
distribution. The Woods-Saxon with A = 5 provides a simi- 
lar case as the Gaussian, both have considerable contributions 
at velocities much higher than the average which will finally 
lead to an inhanced yield at higher transverse momenta for the 
same inverse slope a low pT, i.e. to stronger curvature of the 
spectrum. 

Fig. H shows similar distributions for a high average veloc- 
ity of (Pt) ~ 0.76. Again the box profile leads to a sharp upper 
limit in velocity, but the triangular shape of the distribution is 
distorted by the non-linearity of the relation between Pt and 774 . 
The Woods-Saxon for Z\ = 50 is again similar to the box, and 
the one for A — 5 and the Gaussian are also close. However, 
in this case the two latter distributions have a maximum shifted 
towards higher velocity which is also due to the non-linearity. 
Still, the general property of these distributions is that they have 
larger contributions at high velocities for the same average ve- 
locity as the other two distributions. For such large velocities 
it is noteworthy that the average velocity depends significantly 
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Fig. 2. Velocity distributions as in Fig. |l] at an average velocity of 
{/3t) ~ 0.76. 



on the recipe used to calculate it. If instead of equation ^ one 
uses: 



{PT) = 



(6) 



3 Hadron Spectra from SPS 

The model described above has first been compared to data 
from SPS heavy ion experiments. Data on pion, proton and 
antiproton production in the 3.7 % most central Pb+Pb colli- 
sions from NA44 [|l3|] and in the 10 % most central Pb+Pb 
colUsions from NA49 [ p"4p^ have been used. The trans- 
verse mass spectra of NA44 are shown in Fig. ^ they cover 
the range of O.SGeV/c^ < tot - < 1.3GeV/c^ for pi- 
ous and 0.02GeV/c2 < tot - toq < 0.74GeV/c2 for pro- 
tons and antiprotons. The model has been fitted to the data. 
The best agreement can be achieved with a kinetic tempera- 
ture Tkin — 122.2 McV, an average transverse flow velocity 
(Pj.) = 0.478, a chemical temperature Tchem = 143.9 MeV 
and a baryonic chemical potential fiB = 193.3 MeV. For this 
fit the width parameter has been set to Z\ = 50, allowing for 
a free variation of this parameter produces a slightly smaller 
value with a very large error while the other parameters remain 
unchanged. As the fits seem to be insensitive to small changes 
in A, we have performed most fits with a fixed value. For sys- 
tematic checks some of the constraints on the fits have been 
varied. Most of these changes, as e.g. omitting some of the data 
points (the lowest or highest in rriT, resp.) from the fit, do not 
affect the parameter values significantly. Significant changes 
are obtained under the following conditions: 

1. All contributions from weakly decaying particles are ig- 
nored. 

In this case the chemical temperature and the baryonic 
chemical potential change slightly, while the kinetic tem- 
perature and the flow velocity are essentially unaltered. 
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the numerical values are larger by almost 10%, and differ slightly 
for the different profiles. I will use the definition ^ throughout. 



Fig. 3. Transverse mass distributions of pions and (anti-)protons for 
central Pb+Pb collisions at 158 AGeV from the NA44 experiment 
[p3||. Included are two different fits to the data with the hydrodynamic 
model discussed in the text. 
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Fig. 4. Transverse mass distributions of pions and (anti-)protons for 
central Pb+Pb collisions at 158 AGeV from the NA49 experiment 
[ p4|p^ . Included are two different fits to the data with the hydrody- 
namic model discussed in the text. 



This is understandable as the particle ratios may depend 
more strongly on this contribution than the shape of the 
spectra. As the quality of the fit is similar to the best fit 
above and a not precisely specified fraction of weak decays 
may contribute to the data, this has to be taken into account 
in the systematic error 
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Table 1. Fit parameters of hydrodynamic fits to pion and (anti-)proton spectra in central Pb+Pb collisions at the CERN SPS. The errors given 
are statistical. No error indicates that this parameter was fixed in the fit. 



Remarks 




Tfc™ (MeV) 


■ni 




Tcftem (MeV) 


/is (MeV) 


A 


NA44 data 


best fit 


217/82 


122.2 ± 1.9 


0.805 ± 0.007 


0.478 ± 0.004 


143.9 ±0.6 


193.3 ±0.9 


50 


no weak decays 


217/82 


120.4 ± 1.9 


0.804 ± 0.007 


0.478 ± 0.004 


150.9 ±0.4 


202.6 ±0.9 


50 


broad profile 


360/82 


95.6 ±3.0 


0.745 ±0.010 


0.517 ±0.007 


146.3 ±0.4 


196.4 ±0.9 


5 


Tkin — Tchem 


358/85 


144 


0.728 ± 0.004 


0.442 ± 0.002 


144 


193.3 


50 


high {Pt) 


799/83 


92.7 ±0.9 


0.968 ± 


0.55 


144.6 ± 0.4 


193.2 ±0.9 


50 


high Tchem 


2394/83 


156.3 ±3.0 


0.802 ± 0.001 


0.477 ± 0.001 


165 


208.4 ±0.8 


50 


NA49 data 


best fit 


129/31 


75.8 ±1.9 


0.968 ±0.010 


0.550 ± 0.004 


133.2 ±0.5 


184.5 ±2.0 


50 


low mr omitted 


16/22 


108.0 ±6.6 


0.828 ± 0.034 


0.489 ±0.016 


136.6 ±0.9 


184.6 ±3.3 


50 



2. The width parameter is set to Z\ = 5. 

In this case the chemical parameters show small changes, 
while the temperature and flow velocity change more sig- 
nificantly. This is mainly due to the fact that this spatial dis- 
tribution, which is similar to a Gaussian, results in a broader 
velocity distribution. The quality of the fit is significantly 
worse than the fits given above. 

3. The kinetic and chemical temperatures have been set to the 
same value of T = 144 MeV and the chemical potential to 
/XB = 193.3 MeV. 

Here, the flow velocity obtained is smaller than in the best 
fit, as is expected for a higher kinetic temperature. Again 
the fit quality is worse. 

4. The expansion velocity has been fixed to a value of = 
0.55. The fit quality is much worse than the best fit. 

5. The chemical temperature has been set to Tc/iem = 165 MeV 
similar to results from fits of hadrochemical models to ra- 
tios of total multiplicities of different species [[l6|]. 

This fit cannot describe the spectra, it essentially has the 
wrong ratio of pions to protons. Of importance may be the 
different integration region in rapidity compared to [p^. 
If e.g. the rapidity distribution of protons is broader than 
the one of pions this may easily explain the smaller chemi- 
cal temperature obtained here. However, for the purpose of 
this paper, we will just note this discrepancy and use the 
fit results as they are. The parameters obtained apparently 
describe the particle ratios at midrapidity investigated here, 
and this may be the relevant information which should be 
used also in estimating the contribution of resonances to the 
momentum spectra at this given rapidity. 
Furthermore, it may be that part of the longitudinal mo- 
tion of the source is not generated hydrodynamically but is 
a remnant of the initial state motion, which may be even 
more likely for participant protons. In this case it is not 
completely obvious, whether the protons observed at very 
different rapidities share the same chemical freeze-out tem- 
perature. If the spacetime-momentum coiTelation originates 
from a very early phase of the collision, a unique chemical 
temperature might not be applicable and the use of inte- 
grated multiplicities might be misleading. Rather, there is 
the possibility of a local temperature which may be appli- 
cable for limited rapidity regions. 

The fit parameters obtained from the fits discussed above are 
summarized in table None of the fits provides a perfect agree- 



ment with the data, which is due to a structure in the experimen- 
tal spectra which can not be described within this model. 

The transverse mass spectra of NA49 cover the range 
of 0.05GeV/c2 < tot - mo < 0.75GeV/c2 for pions and 
OGeV/c^ < Tot — toq < 1 GeV/c^ for protons and antipro- 
tons, they are shown in Fig. || together with fits of the model. 
As pions and protons are not available for the same central- 
ity selection - pions have been measured for 5 % central and 
(anti-)protons for 10 % central - the pions have been rescaled 
by a factor of 1/1 -08 according to the number of participants in 
both samples [n7|. This will introduce an additional systematic 
eiTor to the chemical parameters of the fits, which depend most 
significantly on the relative normalization of different species. 
The fitting conditions have been varied as above, the relative 
variation of the parameters is also similarly small. In fact, the 
chemical parameters are relatively similar with slightly smaller 
values compared to the NA44 data. However, for these data the 
best agreement can be achieved with a very low kinetic tem- 
perature Tkin = 75.8 MeV and a very high average transverse 
flow velocity (/3t) = 0.55. It turns out that this significant dif- 
ference to the NA44 data is to a large extent due to the behavior 
of the data (particularly of the protons and antiprotons) at very 
low Tot. If the lowest three data points of all spectra are omit- 
ted from the fit, a kinetic temperature Tkin — 108.0 MeV and 
an average transverse flow velocity {(3t) = 0.489 are obtained 
with a much better x^. The parameters of this fit are relatively 
close to those obtained from the fit to NA44 data. The low tot 
data points of NA49 appear to disagree with the measurements 
by NA44. It is also seen from the best fit to NA49 data that the 
agreement is not perfect if all data points are included - there 
is a tendency of a larger inverse slope of the fit compared to the 
data especially for the antiprotons. 

While the data used above have a limited range in trans- 
verse mass, one might also consider comparing such hydrody- 
namic distributions to data with a wider transverse momentum 
range. This has actually been done before at the SPS by WA98 
[01, using also a program derived from the one of Wiedemann 
and Heinz [^]. However in this version the original assumption 
of the same temperature both for chemical and kinetic freeze- 
out was applied. Also, only neutral pion spectra were used, 
which give very limited information to determine all relevant 
parameters. 

In the following, information from the above fits will be 
used while fitting the neutral pion spectra of WA98 JT^]. Ide- 
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Fig. 5. Transverse mass distributions of neutral pions for the 12.7% 
most central Pb+Pb collisions at 158 AGeY from the WA98 experi- 
ment [ p^ . Included are two different fits to the data with the hydrody- 
namic model discussed in the text. 
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Fig. 6. Transverse mass distributions of neutral pions as in Figure I for 
the 1% most central Pb-l-Pb collisions at 158 AGeV from the WA98 
experiment |[l8]. 



ally one might just try to use exactly the same fits and compare 
it to the neutral pions. As the data are from different experi- 
ments and do not exactly use the same centrality selections, the 
normalization has to be kept as a free parameter The central- 
ity selections used below should however be similar enough, so 
that the shape of the spectra should not differ strongly. 

In Figure ^ fits using the parameters obtained above are 
compared to the neutral pion spectra for the 12.7% most cen- 
tral reactions of Pb-nPb at 158 AGeV. The "best fit" to NA44 in 
table |]is shown as a dashed line, the one to NA49 data as a dot- 
ted line. The description is only superficially adequate. The fit 
to NA44 overshoots the data at large transverse masses, and one 
obtains a x^/^^ = 712/20. The fit to NA49 is below the data 
at large transverse masses with a x^/*^ — 1395/20. This large 
difference for the pion spectra at high rriT is of course due to 
the much larger flow velocity for the NA49 fit, which increases 
the proton inverse slope while decreasing the one of the pions 
compared to NA44. Another fit has been performed keeping 
only the parameters Tchem and /xb the same as in the earlier 
fit and optimizing the kinetic temperature and the flow for the 
neutral pions. This fit (shown as a solid line) yields a good de- 
scription with x^/*^ — 27/18. The fit parameters obtained are 
Tk^n = 130.3 MeV and (Pt) = 0.415. The (dis-)agreement 
of the fits can also be judged from the inset in Figure which 
shows the ratios of the data to the fits. 

One might argue that the different centrality selections are 
responsible for the disagreement. Alt houg h this does not seem 
to be very likely, as it has been shown [|l8|], that the shape of the 
pion spectra does not vary strongly for medium-central to cen- 
tral collisions, a comparison to the 1% most central data has 
been performed. The results are displayed in Figure ^. Again 
the dashed and dotted lines show fits with all parameters fixed 
as above from NA44 and NA49. The behavior of the fits is qual- 



itatively similar to the one discussed above, the disagreement of 
the NA49 fit to the data is not as pronounced as before, which 
is due to the slightly steeper spectra in this case. The solid line 
shows a good fit with Tkin 113.0 McV and (/3t) = 0.459. 
The different parameters for this second fit compared to the 
other central (12.7%) sample are most likely not conclusive, 
but are an example that fitting momentum spectra of just one 
particle species leaves an ambiguity in the parameters. The 
higher flow velocity can be compensated by a lower temper- 
ature and vice versa. Also I will not put too much emphasis on 
the fact that the fits to the NA44 data appear to overpredict the 
pion yield at high px, because part of this discrepancy might 
be due to systematic differences between the two experiments 
involved. 

As one of the major uncertainties in this analysis appears to 
be related to different experimental data sets, it may be worth- 
while to compare data of NA44 and NA49 directly. This has 
been done in Fig. ^by plotting the ratio of the data to the best 
fit of the parameterization to NA44 data (first row in table [l|). 
The same fit is used for both data sets, there is one overall nor- 
malization for all species which differs between experiments. 
One can clearly see that the deviations are most significant at 
low niT for protons and antiprotons. In the systematic stud- 
ies it was found that the parameters for the fits to NA44 do 
not depend very much on the fitting range, while the parame- 
ters for NA49 change significantly when the lowest itit data 
points are not used. It should also be noted that these fits use 
data with statistical eiTors only, as systematic eiTors were not 
readily available for all data sets. If systematic errors would 
account for a significant fraction of the discrepancy between 
the two experiments, a simultaneous fit might be possible and 
yield more conclusive information. This could not be attempted 
in this paper 
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namical parameterization for central Pb+Pb collisions at 158 AGeW 
from NA44 and NA49. 



One can however note that the hydrodynamic parameteri- 
zation can account for a large fraction of the pion production at 
high pt, even if a profile with a negligible tail towards higher 
velocities (Z\ = 50 - similar to the box profile) has been used. 
A more diffuse profile (A = 5) would necessarily lead to a 
larger yield at high px- Most of the other systematic variations 
studied do not lead to significant variations in this respect. The 
most important uncertainty relates to the discrepancy akeady 
seen between NA44 and NA49. Ignoring the very low tot data 
points of NA49 leads to a moderately high flow velocity. In 
this case the pion spectra may be almost entirely be explained 
by such a source. A very high flow velocity as suggested by 
the low Tot NA49 data would considerably reduce the hydro- 
dynamic yield of pions at high tot- 

A fixed chemical temperature of Tchem = 165 MeV leads 
to different results, but from Fig. |3]it is obvious that this would 
also lead to a larger overprediction at high px- As sketched 
above, there may be reasons why this assumption may not be 
adequate for our analysis. In fact, it has been remarked in [ p^pO| ] 
that there is a rapidity dependence of chemical parameters. The 
chemical temperature of Tchem = 143.9 MeV and the baryo- 
chemical potential of hb ~ 193.3 MeV of the best fit can be 
compared to results of other analyses. As stated above, the tem- 
perature is lower than the values of Tchem = 158 ± 3 MeV 
10 and Tchem = 168 ± 2.4 MeV obtained from rapidity 
integrated yields, but cornpares well with the result Tche 
141 ± 5 MeV given in Q. 



4 Hadron Spectra from RHIC 

The situation at RHIC is more favorable for this analysis, as 
there exist data on different particle species partially reaching 
out to large pT measured under the same conditions within one 
experiment. The PHENIX experiment has presented spectra of 
neutral pions [||| and identified charged hadrons | |2^^ in cen- 
tral AuH-Au colHsions at ^snn — 130 GeV. I will also use the 
data from the STAR experiment [|4[^|6[, which only cover 
low Tot. From the STAR experiment, the most recent data on 
protons and antiprotons use a sample of the 6% most central 
collisions, while the pions have been extracted from the 5% 
most central collisions. I will ignore this difference, as it should 
not significantly affect the shapes of the spectra, however, there 
may be small effects on the relative normalization which would 
result in an additional systematic error for the chemical param- 
eters. 

A similar procedure as for the SPS data above is followed, 
with the addition of using pr dependent systematic errors. For 
this a systematic error of 5% added in quadrature to the er- 
rors for the PHENIX data (which contain part of the systematic 
error already) has been assumed [^. For STAR data a sys- 
tematic error of 8% has been used for the antiproton spectra 
[^^. The proton and pion spectra [ ^5p^ were already quoted 
with a systematic error. These does not include all of the nor- 
malization errors, so they might still further affect the chemi- 
cal parameters extracted. Fits are performed to the charged pi- 
ons for Tot — Too < 2 GeV/c and to (anti)protons for tot — 
Too < 3GeV/c. Systematic variations of the fitting assump- 
tions are done as above. In addition I have studied fits to a lim- 
ited transverse mass range tot — too < 1.5 GeV/c for pions 
and (anti)protons in PHENIX. 

The results are summarized in table |[ a comparison of 
some of the fits to the PHENIX data is shown in Fig. ^ those 
to the STAR data are displayed in Fig. |[ 

The fits to PHENIX data show slightly better agreement 
for a wide profile {A = 5), however, the difference in is 
small. As can be seen in Fig. ^ the fit for A = 5 (dotted 
line) has a more important contribution at high pT- The STAR 
data are generally better described using the narrow profile, 
the differences are, however, not easily discernible in Fig. ^ 
For PHENIX, the two different kinds of fits yield very simi- 
lar fit parameters, while there are slightly stronger variations 
when looking at the STAR data. As there does not seem to be 
a strong preference for the wider profile in any of the data sets 
and as there are also theoretical arguments from true hydrody- 
namical calculations [Q, I will base most of the discussion 
below on the fits with a narrow profile. For both data sets fits 
ignoring weak decays yield higher chemical temperatures and 
slightly lower kinetic temperatures. The biggest difference is 
again between the two experiments. The chemical parameters 
still show only a small effect, e.g. while fit (A) to PHENIX data 
yields Tchern 173.3 MeV and = 35.9 MeV, the values 
for STAR (B) are sHghtly smaller with Tchem = 165.3 MeV 
and Us = 31.2 MeV. These values are similar to the results 
in where Tchem = 174 ± 7 MeV and //s = 46 ± 5 MeV 
are given. The better agreement at RHIC compared to SPS of 
the chemical parameters of this analysis with those of the inte- 
grated yields may be related to the fact that boost invariance is 
a much better approximation in the RHIC case. Larger differ- 
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Table 2. Fit parameters of hydrodynamic fits to pion and (anti-)proton spectra in central Au+Au collisions at RHIC. The errors given are fit 
errors using statistical and partial systematic errors (for details see text). No error indicates that this parameter was fixed in the fit. 



Remarks 




Tk^n (MeV) 




(/3t) 


Tchem (MeV) 


MS (MeV) 


A 


PHENIX data 


narrow profile (A) 


101/57 


143.0 ±6.8 


0.745 ± 0.031 


0.450 ±0.016 


173.3 ±2.0 


35.9 ±4.2 


50 


wide profile (A') 


98/57 


142.2 ±7.3 


0.558 ± 0.022 


0.414 ±0.014 


173.0 ±2.0 


37.1 ±4.3 


5 


low rriT - narrow 


80/39 


137.0 ± 7.4 


0.759 ± 0.033 


0.457 ±0.017 


172.5 ±2.0 


35.0 ±4.3 


50 


low rriT - wide 


72/39 


128.5 ±8.3 


0.609 ± 0.026 


0.444 ± 0.015 


172.5 ±2.1 


35.8 ±4.4 


5 


no weak decays - narrow 


116/57 


144.7 ± 7.5 


0.730 ± 0.033 


0.442 ± 0.016 


186.4 ±2.5 


38.6 ±4.6 


50 


no weak decays - wide 


109/57 


141.3 ±8.0 


0.557 ±0.024 


0.413 ±0.015 


185.5 ±2.6 


39.9 ±4.6 


5 


Tkin = Tcheni " narrow 


163/59 


165 


0.633 ± 0.009 


0.392 ± 0.004 


165 


43.0 ±5.1 


50 


Tkin ~ Tchem " widc 


158/59 


165 


0.482 ± 0.006 


0.366 ± 0.004 


165 


43.8 ±5.1 


5 


STAR T - narrow 


186/58 


94.8 


0.962 ± 0.007 


0.548 ± 0.003 


163.6 ± 1.5 


35.1 ±4.1 


50 


STAR T - wide 


206/58 


83.6 


0.730 ± 0.005 


0.509 ± 0.002 


161.4 ± 1.5 


37.5 ±4.1 


5 


STAR data 


narrow profile (B) 


19/36 


94.8 ± 7.9 


1.037 ±0.031 


0.578 ±0.012 


165.3 ±2.2 


31.2 ± 2.5 


50 


wide profile (B') 


36/36 


83.6 ±8.3 


0.984 ± 0.033 


0.622 ±0.013 


174.4 ± 2.9 


33.8 ±2.7 


5 


no weak decays - narrow 


22/36 


88.0 ± 7.1 


1.038 ± 0.030 


0.578 ±0.012 


174.6 ±2.6 


33.0 ±2.7 


50 


no weak decays - wide 


41/36 


77.7 ±7.4 


0.984 ± 0.032 


0.622 ±0.013 


184.9 ±3.3 


36.0 ±2.8 


5 


Tkin = Tchem " narrow 


55/38 


165 


0.920 ± 0.020 


0.530 ± 0.008 


165 


32.4 ±2.5 


50 


Tkin ~ Tchem " widc 


75/38 


165 


0.781 ±0.019 


0.535 ± 0.009 


165 


34.0 ±2.5 


5 


PHENIX T - narrow 


41/37 


143.0 


0.944 ± 0.029 


0.540 ±0.012 


164.4 ±2.4 


32.1 ± 2.5 


50 


PHENIX T - wide 


61/37 


142.2 


0.858 ± 0.032 


0.570 ±0.014 


170.4 ±3.2 


34.2 ± 2.6 


5 


PHENIX + STAR data 


narrow profile (C) 


206/98 


112.3 ±3.3 


0.886 ±0.014 


0.515 ±0.006 


165.0 ± 1.0 


34.5 ± 2.2 


50 


wide profile 


323/98 


108.3 ±4.0 


0.673 ±0.012 


0.479 ± 0.006 


160.9 ± 1.0 


38.1 ± 2.1 


5 


no weak decays - narrow 


230/98 


109.4 ±3.5 


0.888 ±0.015 


0.516 ±0.006 


175.5 ±1.3 


36.7 ±2.3 


50 


no weak decays - wide 


332/98 


103.8 ±4.2 


0.680 ±0.012 


0.483 ± 0.006 


170.4 ± 1.2 


40.3 ± 2.3 


5 


Tkin = T^heni " naiTow 


389/100 


165 


0.672 ± 0.006 


0.413 ±0.003 


165 


38.5 ± 2.1 


50 


Tkin = Tchem " widc 


470/100 


165 


0.517 ±0.004 


0.388 ± 0.002 


165 


40.1 ± 2.1 


5 



ences are seen for the kinetic parameters. Fit (A) yields Tkin = 
143.0 MeV and (/3t) = 0.450, fit (B) Tkm = 94.8 MeV and 
{[3t) = 0.578. As a check also fits have been performed which 
use the kinetic temperature from one experiment as a fixed pa- 
rameter in a fit to the other experiment, and vice versa. A rea- 
sonable fit to the STAR data (with a higher y^) can be obtained 
this way by using a high kinetic temperature, but the flow ve- 
locity turns out to be larger than for PHENIX. In the other 
case, the fit to PHENIX data with a low temperature yields a 
much higher - no good fit can be obtained. There have been 
fits to the PHENIX data using hydrodynamic parameterizations 
that do not include resonances [|^. These yield similar aver- 
age flow {{Pt) — 0.47 ± 0.01) velocities but smaller kinetic 
freeze-out temperatures {Tkin — 121 ± 4 MeV) - this differ- 
ence is consistent with the influence of resonance decays. 

Fits assuming a unique temperature of T = 165 MeV for 
both kinetic and chemical freeze-out as in are again much 
worse in fit quality for both shapes of the source distribution, 
as can be seen from the large in table |[ The assumption of 
simultaneous kinetic and chemical freeze-out underlying this 
fit is questionable, but I will still consider such a fit once more 
below. One might argue that the transverse mass range used in 
these fits does include possible contributions of hard scattering. 
However, fits which use data only for low rriT {rriT — mo < 
1.5 GeV/c) yield very similar results as seen in table ^ 

Finally, a simultaneous fit to STAR and PHENIX data has 
been performed, where the STAR data have been rescaled by 



a factor 0.92 (ratio of the number of participants) to account 
for the different centrality selection. Please note that such a si- 
multaneous fit is only successful if one includes the system- 
atic errors. Considering just statistical errors, the data from 
the two experiments could not be consistently described by 
this model. The results are also given in table ||, the best fit 
(C) is compared to the data in Fig. |l^ Fit (C) is obtained 
using Z\ = 50 and describes the data reasonably well with 
X^/i^ — 206/98. Chemical parameters change relatively lit- 
tle. As might be expected the kinetic parameters are interme- 
diate between those for the individual experiments, one finds 
Tfcm = 112.3 MeV and (/3t) = 0.515. Once again a fit using 
Tkin — Tchem = 165 MeV is significantly worse. 



5 Discussion 

The results for SPS and RHIC data may now be compared 
to each other. In both cases a box-like profile is favored over 
a Gaussian-like profile. It was akeady discussed in that a 
box-like profile would be required in case of a large flow ve- 
locity, and this was also found to be the preferred profile in [|^]. 
Such a distribution with little or no "tail" to higher velocities 
is also closer to expectations from true hydrodynamic calcula- 
tions [^. In both cases (SPS and RHIC) there are fits with a 
range of flow velocities, depending mainly on the experiment 
that one concentrates on. At the SPS the range of velocities is 
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Au+Au 
b-130 AGeV 
1 % central 
(PHENIX) 



12 3 
m^-mp (GeV/c2) 



2 3 4 
m^-mg (GeV/c2) 



Fig. 8. Transverse mass distributions of pions and (anti)protons for the 
10% most central Au+Au collisions at ^/sNN = 130 GeV from the 
PHENIX experiment [|j2|,|3||. The lines show fits (A) and (A') from 
table ^. The inset shows the ratio of the data to fit (A). 
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Fig. 9. Transverse mass distributions of protons and antiprotons for 
the 6% most central Au+Au collisions at ^/sNN = 130 GeV and 
for negative pions for the 5% most central Au+Au collisions from the 
STAR experiment The lines show fits (B) and (B') from 

table p| The inset shows the ratio of the data to fit (B). 



(Pt) — 0.48 — 0.55 and the analysis presented here does not 
allow to constrain this value further. At RHIC a combined fit is 
possible, which yields as the most likely value {(3t) — 0.515. 
Similarly the range of kinetic freeze-out temperatures for the 
SPS is Tkin = 76 - 122 MeV, while the RHIC value appears 
to be Tkin = 112.3 MeV. Apparently this analysis indicates 
that the flow velocities at SPS and RHIC may be similar, and 
there is no clear indication of stronger collective flow at RHIC 
as found in ||3^. It can, however, be seen in that this is 
strongly influenced by the proton and antiproton spectra from 
STAR used in the analysis, which show a much larger inverse 
slope than the PHENIX results. 

Both for SPS and RHIC there are still considerable uncer- 
tainties which are related to the experimental data themselves. 
For the SPS these imply large uncertainties for the extrapola- 
tion of the hydrodynamical fits to higher my- Still, if we at- 
tempt to ascribe particle production at low niT in heavy-ion 
collisions at the SPS to hydrodynamic scenarios, the extrapola- 
tion of the hydrodynamic contribution to high tot will yield a 
considerable contribution of hydrodynamic production. 

For RHIC the measurement of protons and antiprotons at 
higher rriT does provide additional constraints and excludes 
some of the extreme scenarios. A simultaneous fit to STAR and 
PHENIX has been performed, which yields a reasonable agree- 
ment with the available data. In the following I will assume that 
this fit describes the contribution of a hydrodynamic source to 
particle emission. It does describe the spectra almost entirely 
at low niT- At higher rriT the fit starts to deviate, which is not 
surprising, as one would expect additional production mecha- 
nisms, especially hard scattering, to play a role there. However, 
this hydrodynamic contribution may be still significant even at 
high momenta - this would be of particular relevance to any 



study of hard scattering, because to be able to obtain informa- 
tion on hard scattering production one would have to know this 
hydrodynamic contribution and take it into account. 

A simple partitioning of the produced particles into either 
hydrodynamic or hard scattering production is not possible. If 
one believes that equilibration is taking place in these colli- 
sions, there will almost necessarily be parts of the reaction vol- 
ume that are only partially equilibrated. At RHIC the high pT 
part of the non-equilibrated spectrum (with a power-law shape 
from hard scattering) should have a larger inverse slope than 
the equilibrated spectrum (similar to exponential). So in this 
case equilibration through multiple collisions will slow parti- 
cles down and soften the spectrum. As the degree of equili- 
bration will vary throughout the reaction volume, the spectrum 
might be a superposition of a range of spectra between the one 
from the initial collisions and a thermalized one. 

At this point one does not necessarily need to distinguish 
between parton and hadron degrees of freedom. Thus, parton 
energy loss would just be one possible equilibration mecha- 
nism. While it is not possible to precisely predict how close the 
final distribution would be to a thermal (hydrodynamic) distri- 
bution, the latter can be taken as a limiting case. If hard scat- 
tered partons are completely quenched, one will still end up 
with a thermal distribution as a minimum yield. Any excess 
over this minimum might be due to non-equilibrated contribu- 
tions, either true hard scattering production or not completely 
equilibrated production, which still has a higher average mo- 
mentum than the thermal. 

The discussion in this paper assumes that the spectra in the 
intermediate momentum region can be described by thermal 
distributions. This is not completely correct, as e.g. even hard 
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12 3 

mj-mg (GeV/c^) 



2 3 4 
rn^-mg (GeV/c2) 



Fig. 10. Transverse mass distributions of pions and (anti)protons for 
central Au+Au collisions at ^sjvjv = 130 GeV from the PHENIX 
and STAR experiments. The lines show a simultaneous fit (C) from 
table ^. The inset shows the ratio of the data to the fit. 




p-^ (GeV/c^ 



Fig. 11. Fraction of the hydrodynamic component of the neutral pion 
and proton+antiproton yield using fit C from table ^. The dotted line 
shows a fit of a Woods-Saxon function to the data (see text). 



scatterings will produce some lower momentum component. 
However, as any comparison of pure hydrodynamics to data, I 
will assume for the following that this non-thermal contribution 
at low and intermediate pT is negligible. 

In Fig. |ri| ratios of the simultaneous fit (C) to the exper- 
imental data, i.e. the fractions Rhydro of the produced parti- 
cles accounted for by the hydrodynamic parameterization, are 



shown for neutral pions (solid squares) and for protonsH-antiprotons 
(open circles) as a function of px- 

The protonH-antiproton yield is reasonably described by this 
fit in the entire pT range investigated, as can be seen by the 
more or less constant ratio. The hydrodynamic contribution to 
pions depends more strongly on pT- While here hydrodynam- 
ics dominates also at lower pT, there is a steady decrease of the 
contribution towards higher pT- Still even for pt > 3 GeV/c 
the hydrodynamic contribution to the pion yield is close to 
40%. 

The not completely thermalized part of the emitted hadrons 
will reflect itself in other observables. A pT dependent non- 
thermalized contribution to the spectrum will be related to a 
Pt dependent deviation of the measured elliptic flow V2 from 
the purely hydrodynamic behavior, because the latter relies on 
equilibration achieved earlier in the collision history than ki- 
netic freeze-out. Quantitative estimates of this are difficult, as 
the partially thermalized contribution will lead to an unknown 
contribution to V2- Full hydrodynamic calculations, when ap- 
plied to elliptic flow and particle spectra, should be able to ex- 
tract more information on this. 



6 Summary 

Results of hydrodynamic fits to momentum spectra at mid- 
rapidity in heavy ion reactions at the SPS and at RHIC have 
been presented. While at RHIC a simultaneous description of 
STAR and PHENIX data is possible, when systematic errors 
are taken into account, there are still unresolved differences 
at SPS depending on the experiment used. At SPS a chemi- 
cal temperature of Tchem ~ 145 MeV and a baryo-chemical 
potential of /is w 190 McV have been obtained. This chem- 
ical temperature is slightly smaller than values obtained from 
rapidity integrated yields. At RHIC the corresponding values 
are Tchem ~ 165 MeV and ~ 35 McV. The kinetic param- 
eters are Tkm = 76 - 122 McV and {Pt) = 0.48 - 0.55 at 
SPS and Tfc™ w 110 McV and (/3t> ~ 0.52 at RHIC when us- 
ing a velocity profile similar to a box profile. The momentum 
spectra at RHIC are broader than at the SPS. However, this 
analysis does not prove a significantly stronger transverse flow 
at RHIC, the change in the momentum spectra might also be 
due to a higher freeze-out temperature. At high transverse mo- 
menta the spectra of protons can be nicely described at RHIC. 
The pion spectra leave room for additional particle production 
mechanisms (e.g. hard scattering) at high px- Estimates of the 
suppression of hard scattering should, however, take the non- 
negligible contribution from a hydrodynamic source into ac- 
count. Important hydrodynamic contributions in spectra at high 
Pt require even stronger quenching of hard scattering produc- 
tion. In return, quenching is a necessary condition to achieve 
locally equilibrated distributions. 



The author would like to thank U.A. Wiedemann for helpful dis- 
cussions and J.M. Burward-Hoy, M. Calderon de la Barca Sanchez, 
M. Gazdzicki, M. Kaneta, K. Schweda, R. Snellings for valuable in- 
formation regarding the experimental data. Part of this work was cre- 
ated at the University of Miinster and at GSI Darmstadt. 



10 



T. Peitzmann: Hydrodynamic description of spectra at iiigli transverse mass in ultrarelativistic heavy ion collisions 



References 

1. E. Schnedermann, J. Sollfrank, and U. Heinz, Phys. Rev. C 48 
(1993) 2462. 

2. S. Ciiapman, J. R. Nix and U. Heinz, Vhys. Rev. C 52 (1995) 
2694. 

3. P .Braun-Munzinger et al., Piiys. Lett. B 344 (1995) 43. 

4. U.A. Wiedemann and U. Heinz, Phys. Rev. C 56 (1997) 3265. 

5. NA49 Collaboration, H. Appelshauser et al., Eur. Phys. J. C 2 

(1998) 661-670. 

6. WA80 Collaboration, R. Albrecht et al., Eur. Phys. J. C 5 (1998) 
255. 

7. WA98 Collaboration, M.M. Aggarwal et al., Phys. Rev. Lett 83 

(1999) 926. 

8. B. Tomasi k. U.A. Wiedemann, and U. Heinz, preprint n ucl- 
th/9907096. 

9. PHENIX Collaboration, K. Adcox et al., Phys. Rev. Lett. 88 
(2002) 022301. 

10. STAR Collaboration, C. Adler et al., Phys. Rev Lett. 87 (2001) 
182301. 

11. PHENIX Coll aboration, K. Adc ox et al., submitted to Phys. Rev 



Lett., preprint 



nucl-ex/0204005 



12. NA49 Collaboration, RG. Jones et al., Nucl. Phys. A 610 (1996) 
188c. 

13. NA44 Coll aboration, l.G. Be arden et al., submitted to Phys. Rev. 



C, preprint nucl-ex/0205002 



14. NA49 Collabor ation, S.V. Afana siev et al., submitted to Phys 



Rev. C, preprint nucl-ex/0202019 



15. NA49 Collaboration, S.V. Afanasiev et al., proceedings of Quark 
Matter 2002, Nantes, France. 

16. P. Braun-Munzinger, I. Heppe, and J. Stachel, Phys. Lett. B 465 
(1999) 15-20. 

17. M. Gazdzicki, private communication. 

18. WA98 Collaboration, M.M. Aggarwal et al., Eur Phys. J. C 23 
(2002) 225-236. 

19. J. Sollfrank, U. Heinz, H. Sorge, and N. Xu, Phys. Rev. C 59 
(1999) 1637. 

20. M. Kaneta, and N. Xu, J. Phys. G 27 (2001) 589. 

21. F. Becattini, J. Cleymans, A. Keranen, E. Suhonen, and 
K. Redlich, Phys. Rev. C 64 (2001) 024901. 

22. PHENIX Collaboration, K. Adcox et al., Phys. Rev Lett. 88 
(2002) 242301. 

23. PHENIX Collaboration, J. Velkovska et al., talk at ICPA-QGP 
2001 Jaipur, to be published in Pramana - Journal of Physics. 

24. STAR Collaboration, C. Adler et al., Phys. Rev. Lett. 87 (2001) 
262302. 

25. STAR Collaboration, K. Schweda et al., proceedings of 
RIKEN BNL Workshop Vol.42 (2002) p.57,eds. M. Gyulassy, 
D. Kharzeev, and N. Xu, Brookhaven National Laboratory, Up- 
ton, NY 11973-500, USA. 

26. M. Calderonj_^hD_^hesiSj_^ale_^niversi^^ 
ex/0111004. 

27. PHENIX Collaboration, K. Adcox et al., in preparation. 

28. P. Huovinen, P.P. Kolb, U. Heinz, P.V. Ruuskanen, and 
S.A. Voloshin, Phys. Lett. B 503 (2001) 58-64 

29. P. Braun-Munzinger, D. Magestro, K. Redlich and J. Stachel, 
Phys. Lett. B 518 (2001) 41-46. 

30. PHENIX collaboration, J.M. Burward-Hoy et al., talk at Quark 
Matter 2002, Nan tes, to be published in Nucl. Phys. A, preprint 



nucl-ex/02 10001 



31. W. Broniowski and W. Florkowski, Phys.Rev.Lett. 87 (2001) 
272302. 

32. N. Xu and M. Kaneta, Nucl. Phys. A 698 (2002) 306c. 



